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Abstract 

The dynamic impedance of a sphere oscillating in an elastic medium is considered. Oestreicher's pQ 
formula for the impedance of a sphere bonded to the surrounding medium can be expressed simply in terms 
of three lumped impedances associated with the displaced mass and the longitudinal and transverse waves. 
If the surface of the sphere slips while the normal velocity remains continuous, the impedance formula is 
modified by adjusting the definition of the transverse impedance to include the interfacial impedance. 
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The dynamic impedance of a spherical particle embedded in a medium is important for acoustical measure- 
ment and imaging. The impedance is used, for instance, in measurement of the mechanical properties of 
tissue, e.g. [2], and is intimately related to the radiation forcing on particles The latter is the basis 
for imaging techniques such as vibro-acoustography which has considerable potential in mammography for 
detection of microcalcifications in breast tissue ^ . Oestreicher derived the impedance for a rigid sphere 
oscillating in a viscoclastic medium over 50 years ago. Although it was derived for a sphere in an infinite 
medium, Oestreicher's formula is also applicable, with minor modification j5|, to dynamical indentation 
techniques where the particle is in contact with the surface of the specimen, see Zhang et al. [H] for a review 
of related work. Chen et al. recently validated Oestreicher's formula experimentally by measuring the 
dynamic radiation force on a sphere in a fluid. The impedance formula is based on perfect no slip conditions 
between the spherical inclusion and its surroundings. This may not always be a valid assumption, e.g., in 
circumstances where a foreign object is embedded in soft material. This was precisely the situation in recent 
measurements of the viscosity of DNA cross- linked gels by magnetic forcing on a small steel sphere jHj . This 
paper generalizes the impedance formula to include the possibility of dynamic slip. 

Two related results are derived in this paper. The first is a modified form of Oestreicher's formula which 
enables it to be interpreted in terms of lumped parameter impedances. This leads to a simple means to 
consider the more general case of a sphere oscillating in a viscoclastic medium which is permitted to slip 
relative to its surroundings. The slip is characterized by an interfacial impedance which relates the shear 
stress to the discontinuity in tangential velocities. This generalization includes Oestreicher's original formula 
as the limit of infinite interfacial impedance, and agrees with previous results for the static stiffness of a 
spherical inclusion with and without slip 0. 

2 Summary of results 

A sphere undergoes time harmonic oscillatory motion of amplitude uq in the direction x, 

u sphere = u e- iut ± (1) 

The time harmonic factor e~ lu)t is omitted but understood in future expressions. The sphere, which is 
assumed to be rigid and of radius a, is embedded in an elastic medium of infinite extent with mass density p 
and Lame moduli A and (i. The moduli may be real or complex, corresponding to an elastic or viscoelastic 
solid. We will later consider complex shear modulus fi = \i\ — 110^2, where the imaginary term dominates in 
a viscous medium. The force exerted on the sphere by the surrounding medium acts in the x-direction, and 
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is defined by 



where T is the traction vector on the surface. The sphere impedance is defined 

F 



Fx= / Tds, 
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Oestreicher's expression for the impedance of a sphere that does not slip relative to the elastic medium is 1 
□ 
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Here k and /i are, respectively, the longitudinal and transverse wavenumbers, k = u/cl, h = lo/ct with 
c L = \/(X + 2fi)/p and c T = \fp~Jp. 

Noting that Oestreicher's formula can be rewritten 

i -i - 



Z = —na piu) ^ — 1 + 
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implics our first result, that the impedance satisfies 
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Z + Z m Zl + Z m Zt + Z m ' 
where the three additional impedances are defined as 

7 ■ 4 3 

Z L = (iw) _1 47ra(A + 2p) (1 - ika), 
Zt = (iuj)^ 1 47ra/i(l — iha). 



(5) 
(6) 

(7a) 
(7b) 
(7c) 



The second result is that if the sphere is allowed to slip relative to the elastic medium then the general 
form of eq. © is preserved with Zt modified. Specifically, suppose the tangential component of the traction 
satisfies 



T • t = zi (u s f /lere - v) • t, 



(8) 



where t is a unit tangent vector, v the velocity of the elastic medium adjacent to the sphere, and zj is an 
interfacial impedance 2 , discussed later. Equation JSJ) holds at each point on the interface r = a. We find 
that Z now satisfies 

1 1 9 

(9) 



Z + Z n 



Zl + Z m Zs + Z n 



1 Equation J3} is Oestreicher's eq. (18) with i replaced by — t since he used time dependence e iu *. 

2 Capital Z and lower case z are used to distinguish impedances defined by force and stress, respectively. 
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where the new impedance Z$ is given by 

1 1 1 



(10) 



Zs Zt 47ra 2 z/ + (iu) 1 87ra/z 
These results are derived in the next Section and discussed in Section 0] 

3 Analysis 

We use Oestreicher's P representation for the elastic field outside the sphere, 

u = -A x grad ( x ) + Bi[2h {hr) gradz - h 2 (hr)r 3 grad -^] , r > a, (11) 

where r = |r| is the spherical radius and x is the component of r in the x— direction, both with origin at the 
center of the sphere. Also, h n are spherical Hankel functions of the first kind ^U|- Let f = r~ 1 r denote the 
unit radial vector, then 

u = -Ax [MM X _ h 2 (kr)-r] + Bi \2h (hr)St - h 2 (hr){± - 3-r)l . (12) 
kr r r 

The surface traction is T = erf where er is the stress tensor. The traction can be calculated from Ijl2|l and 

the following identity pQ for an isotropic solid, 



l d 1 

T = f Adiv u + — gradr • u + )u. (13) 

r or r 



Thus, referring to @, we have 



2 2 a 

T •x=[2/i/i 2 (fcr)(l- 3^) +(A + 2/j,)krh 1 (kr) ^] -± 

+ [2h 2 (hr)(l - 3^J) - hrhx(hr)(l - ^)] 3 M ^ . (14) 

Integrating over the sphere surface, the resultant is 

„ 4 , or, h-i(ka) „ hi(ha)-, 

F=-na 3 pu ) 2 [A 1 ^^-6B 1 ^j—L}. (15) 

The coefficients A\ and B\ follow from the conditions describing the interaction of the sphere with its 
surroundings. These are the general slip condition (JSJ) plus the requirement that the normal velocity is 
continuous. The conditions at the surface of the sphere are 

u ■ r = uq x • f 
T • t = iuzi (u — uox) • t 

By symmetry, the only non-zero tangential component is in the plane of f and x, and we therefore set 
t = 8 = (f cos# — x)/ sin^ where = arccosf ■ x is the spherical polar angle. Using polar coordinates, 
u = u r r + ugO and T = a rr r + <r r e9, and (|16fl becomes 

u r = uq cos 9 I 
<r r g — iuizjug = iu!Zju sin J 



r = a. (16) 



r = a, 0<6 <n. (17) 
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The shear stress follows from the identity 



,du e ldu r , 1oN 
or r oO r ' 

and the interface conditions (|17|) then imply, respectively, 

[^a(fco) - ^] Ai + 6 ^ ' B 1= uo. (19a) 

r h-\(ka) 2rih2(ka) 1 , r , uha\hi(ha) . 2a , n „ , 

- + Ai + (2 + ? )^t— L - 1 + — — )/»2 H 3fli = uo. 19b 

feet luazj lujazj na itoazj 

Solving for Ai and B\, then substituting them into eqs. (|15|l and @, and using the known forms for the 
spherical Hankel functions, yields 

Z = -Z m + 3/[l/(Zi + Z m ) + 2/(Z s + Z m )] . (20) 

Equation l|20|) is identical to I©, which completes the derivation of the generalized impedance formula. 

4 Discussion 

It is useful to recall some basic properties of lumped parameter impedances. The impedance of a spring 
mass damper system of stiffness K, mass M and damping C is 

Z= (ito^K -C + iwM. (21) 

Two impedances Z\ and Z^ combined in series have an effective impedance (Z^ 1 + Z^ 1 )^ 1 , while the result 
for the same pair in parallel is [Z\ + Z2). 

Referring to the definitions of cq. I|7)l. it is clear that Z m is the impedance of the mass of the volume 
removed from the elastic medium. The impedance of a longitudinal or transverse plane wave is defined as 
the ratio of the stress (normal or shear) to particle velocity, and equals zl, Zt, where 

zl = PCl, zt = per- (22) 

Thus, both Zl and Zt have the form 

Z= (— - l)47ra 2 z, (23) 
ina 

where k is the wavenumber (k or h). In particular, the impedances Zl and Zt have stiffness and damping, 
but no mass contribution. The damping can be ascribed to the radiation of longitudinal and transverse 
waves from the sphere. 

The impedance Z§ of eq. (|l(Jf) corresponds to Zt in series with an impedance Zj, where 

Zi = 47ra 2 2 7 + — - . (24) 

iha 
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Thus, Zi can be interpreted as the total interfacial impedance for the surface area of the sphere in parallel 
with twice the stiffness part of Z?- 

The limit of a purely acoustic fluid is obtained by letting the shear modulus (i tend to zero with A finite, 
while an incompressible elastic or viscous medium is obtained in the limit as the bulk modulus A + |/x 
becomes infinite with [i finite. The acoustic and incompressible limits follow from H20|) as 

acoustic medium, 

(25) 

i \Zg 1 incompressible medium. 

Thus, Z for the acoustic fluid comprises \Z Ta in series with \Z^. Note that, as expected, the interfacial 
impedance zi is redundant in the acoustic limit. The impedance for the incompressible medium is \Z m and 
\Zs in parallel, and it depends upon the interfacial impedance. 

In order to examine the role of zi , we first express the impedance Z of eq. i|2(J|) in a form similar to J5J , 




Z ~ —7ra 3 piuj I — 1 
3 H 



1, 3(1 — ika)^ -i 2 , 3(1 — iha) 



3^ k 2 a 2 ' ' 3 V ~ h 2 a 2 [l + ^{l-iha)]' 



(26) 



2 ' 

where the influence of the interfacial impedance is represented through the non-dimensional parameter 

X=(l + =)" 1 - (27) 

The form of \ is chosen so that it takes on the values zero or unity in the limit that the sphere is perfectly 
bonded or is perfectly lubricated, 

fo, no slip, *,^oo, 
I 1, slip, zi = 0. 

The acoustic and incompressible limits of l|25(l are explicitly 

3 ™ P IUJ 2iFitoFW . acoustic, 

(29) 

^ [ i+W^ha) ~ ] - incompressible. 

Oestreicher £Q showed that the original formula Q provides the acoustic and incompressible limits for perfect 
bonding (x = 0). Ilinskii ct al. |1 1 j derived the impedance in the context of incompressible elasticity, also 
for the case of no slip. 

The behavior of Z at low and high frequencies depends upon how zj and hence x behaves in these limits. 
For simplicity, let us consider \ as constant in each limit, equal to xo at low frequency, and Xoo at high 
frequency. Then, 

Z = <! (30) 
lna 2 pc L [ - (1 + 2£(1 - Xoo)) + ^(l - 4£[1 + - l) Xoo ]) 

+o( F ^)], |H,M»i- 



Impedance of an oscillating sphere February 6, 2008 7 

The leading order term at high frequency is a damping, associated with radiation from the sphere. The 
dominant effect at low frequency is, as one might expect, a stiffness, with the second term a damping. The 
low frequency stiffness is identical to that previously determined by Lin et al [§] who considered the static 
problem of a sphere in an elastic medium with an applied force. They derived the resulting displacement, 
and hence stiffness, under slip and no slip conditions. In order to compare with their results, we rewrite the 
leading order term as 

5 — w + 2(1 - v)xq 



where v is the Poisson's ratio, 



2°L 



(32) 

Equation i|31|) with \o = and xo = 1 agrees with eqs. (40) and (41) of Lin et al. [5], respectively. In an 
incompressible viscous medium with v « 1/2 and \x = —iu>fi2, <|3 1 1) becomes 

(33) 

1 + -5X0 

which reduces to the Stokes drag formula F = —Qira^v for perfect bonding. When there is slip (\o = 1) 
the drag is reduced by one third, F = —Aira^v- It is interesting to note that one third of the contribution to 
the drag in Stokes' formula is from pressure, lim^v^ the remained from shear acting on the sphere. However, 
under slip conditions, the shear force is absent and the total drag Ana^v is caused by the pressure. 

The simplest example of the interfacial impedance is a constant value, which is necessarily negative and 
corresponds to a damping, zj = — C . For an elastic medium we have 

X = ; — , w c = — j-, elastic medium, zj = —C. (34) 

1 — iuj/uj c aC 

Hence, xo = 1 and %oo = 0, corresponding to slip at low frequency and no slip at high frequency. The 
transition from the low to high frequency regime occurs for frequencies in the range of a characteristic 
frequency uj c . Alternatively, if the medium is purely viscous fi = —iu>H2, again with constant zj, the 
parameter x becomes 

X = aC~i viscous medium, /i = —icjfj,2, zj = —C. (35) 

1 + 2717 

In this case x is constant with a value between and 1 that depends upon the ratio of the interfacial to bulk 
viscous damping coefficients, and also upon a. One can define a characteristic particle size a c = ^/C ', such 
that spheres of radius a <C ac ("> ac) arG effectively bonded (lubricated). 

Figures 1 and 2 show the reactance and resistance of a sphere of radius 0.01 m in a medium with 
the parameters considered by Ocstrcicher [P based on measurements of human tissue, p = 1100 kg/m 3 , 
/ii = 2.5 x 10 3 Pa, (12 = 15 Pa sec, A = 2.6 x 10 9 Pa. The perfectly bonded (x = 0) and perfect slip (x = 1) 
conditions are compared. Figure 1 indicates that the mass-like reactance is generally reduced by the slipping, 
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and it also shows that the low frequency stiffness is two-thirds that of the bonded case, eq. (|31|l . Interfacial 
slip leads to a significant decrease in the resistance, as evident from Figure 2 which shows a reduction for all 
frequencies. 
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Frequency in c.p.s. 



Figure 1: The reactance (real part of Z) for an oscillating sphere in a tissue- like material PQ. The solid and 
dashed curves correspond to a bonded (x = 0) and slipping (x = 1) spherical interface, respectively. The 
reactance is positive (mass-like) except for frequency below 30 c.p.s. (50 c.p.s. for the dashed curve) where 
it is negative (stiffness- like). 
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Figure 2: The resistance (imaginary part of —Z) for an oscillating sphere in a tissue-like material p^. The 
solid and dashed curves correspond to a bonded (x = 0) and slipping (x = 1) spherical interface. 



